Collisional effects on the collective laser cooling of trapped bosonic gases 
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We analyse the effects of atom-atom collisions on collective 
laser cooling scheme. We derive a quantum Master equation 
which describes the laser cooling in presence of atom-atom 
collisions in the weak-condensation regime. Using such equa- 
tion, we perform Monte Carlo simulations of the population 
dynamics in one and three dimensions. We observe that the 
ground-state laser-induced condensation is maintained in the 
presence of collisions. Laser cooling causes a transition from 
a Bose-Einstein distribution describing coUisionally induced 
equilibrium,to a distribution with an effective zero tempera- 
ture. We analyse also the effects of atom-atom collisions on 
the cooling into an excited state of the trap. 



I. INTRODUCTION 

In the recent years, laser cooling has constituted one 
of the most active research fields in atomic physics [nl. 
However, the laser cooling techniques by themselves have 
not allowed to reach temperatures for which the quan- 
tum statistical effects become evident. In particular, only 
the combination of laser cooling, and evaporative M or 
sympathetic cooling Q has permited in the last years to 
observe experimentally, the Bose-Einstein condensation 
(BEC) in alkali gases, seventy years after its theoretical 
prediction B. The question whether it is or it is not 
possible to achieve the BEC only with laser cooling tech- 
niques remain, at least as an intellectual challenge. The 
laser-induced BEC is, however, not only an academic 
problem, but has several advantages with respect to the 
nowadays widely-employed collisional mechanisms (as 
evaporative cooling). These advantages are: (i) the num- 
ber of atoms does not decrease during the cooling pro- 
cess; (ii) it is possible to design a non-destructive BEC 
detection by fluorescence measurements; (iii) reacher ef- 
fects can appear, since now the system is open, i.e. it 
is not in thermal equilibrium; (iv) laser-induced conden- 
sation can be used to design techniques to pump atoms 
into the condensate. The latter is specially important in 
the contex of future atom-laser devices ^-g| . 

The main problem which prevents experimentalists 
from obtaining BEC by optical means is the reabsorp- 
tion of spontaneously emitted photons. The most ef- 
fective laser-cooling techniques (such as VSCPT [|j, or 
Raman cooling jl^), are based on the crucial concept of 
dark states O], i.e. states which cannot absorb the laser 
light, but can receive population via incoherent pump- 
ing, i.e. via spontaneous emission. Unfortunately, the 
atoms occupying the dark states are not unaffected by 



the photons spontaneously emitted by other atoms. This 
problem turns to be very important at high densities, 
as those required for the BEC [Q; in such conditions 
dark-state cooling techniques cease to work adequately. 
Several remedies to the reabsorption problem have been 
proposed, as the reduction of the dimensionality of the 
trap from three to two or one dimensions |Q , or the use 
of traps with frequencies, w, of the order of the recoil fre- 
quency {ujR = 'hk\/2M , where fc^ is the laser wavevector 
and M is the atomic mass) [|6|. Other, perhaps more 
promising, idea consists in exploiting the dependence of 
the reabsorption probability on the fluorescence rate 7. 
In particular, in the so-called Festina Lente limit p4[ , 
when 7 < cj with uj the trap frequency, the heating effects 
of the reabsorption can be neglected. Another proposal 
consists in working in the so-called Bosonic Accumula- 
tion Regime S, in which the reabsorption can, under 
certain conditions, even help to build up the condensate. 
In the following we shall assume that the considered sys- 
tem fulfills the Festina Lente limit. 

In a series of papers ||l^,|lg , we have proposed a cool- 
ing mechanism (which we have called Dynamical cooling) 
which permits the cooling of an atomic sample into an 
arbitrary single state of an harmonic trap, beyond the 
Lamb-Dicke limit (i.e. when the Lamb-Dicke parameter 
1] > 1, with T]^ = lor/llj). The cooling mechanism em- 
ploys laser pulses of different frequencies (and eventually 
different directions, phases and intensities), in such a way 
that a particular state of the trap remains dark during 
the cooling process, acting as a trapping state. There- 
fore, the population is finally transferred to this particu- 
lar state. We have first analysed the particular situation 
of a single atom in the trap 1 15|, and extended the analy- 
sis to a collection of trapped bosons [1^ . We have shown 
that the bosonic statistics helps to achieve more robust 
and rapid condensation, as well as to produce non-linear 
effects, such as hysteresis and multistability phenomena. 

However, all the calculations performed so far in the 
analysis of the dynamical cooling scheme do not take into 
account the atom-atom collisions, i.e. are considered in 
the so-called ideal gas limit. The ideal gas limit imposes 
important restrictions to the physical system, in partic- 
ular the atomic density cannot be very large. An in- 
teresting possibility in order to achieve quasi-ideal gases 
consists in the "switching-off" of the s-wave scattering 
length a (which is the main contribution to the atom- 
atom collisions for sufficiently low energies), either by 
employing magnetic fields (tuning the so-called Feshbach 
resonances |17|] ) or by using a red-detuned laser tuned 
between molecular resonances as proposed by Fedichev 
et al |18| . However, without special precautions, the ef- 



fects of the atom-atom collisions play a substantial role. 
It is the aim of this paper to analyse such effects in the 
context of our dynamical laser cooling scheme. 

In recent years, C. W. Gardiner, P. ZoUer and collab- 
orators have devoted a series of papers [19 - 23 1 to the de- 
cription of interacting Bose gases with and without trap- 
ping potentials. These authors have developed a quan- 
tum kinetic theory of Bose gases. In particular, for the 
case of a weakly interacting gas, a so-called Quantum 
Kinetic Master Equation (QKME) has been formulated 
p9[ , which is a quantum stochastic equation for the ki- 
netics of the dilute Bose gas, that describes the behavior 
and formation of the condensate. This equation is very 
difficult to simulate, and therefore various simplifications 
have been proposed. Particularly interesting results are 
obtaining by using the so-called Quantum Boltzmann 
Master Equation (QBME) which neglects all spatial in- 
homogeinity of the trapped states ||2^ . Although this is 
an extreme simplification, the (very much easier) simula- 
tion of the QBME give a good idea of the solutions that 
the QKME could produce. In the following, we shall show 
that the master equation (ME) which describes the laser 
cooling problem in the presence of atom-atom collisions 
can be, in the case of the weak-condensation regime, 
splitted into two independent parts, one accounting for 
the collisional effects (which has the form of the QBME 
proposed in Ref. 1131), and another which describes the 
laser cooling process, and has the form of the ME already 
developed for the case without collisions |I6[ . 

The structure of the paper is as follows. In Sec. ||, we 
derive the quantum ME which describes the laser cooling 
plus collisions in the weak-condensation regime. In Sec. 
Ill, we present the results for one-dimensional excited- 
state cooling. Sec. |^is devoted to the three-dimensional 
results for the case of ground-state cooling. Here, we 
use additional ergodic approximation which assures fast 
redistribution of atoms within an energy shell. Finally, 
in Sec. M we summarize some conclusions. 



II. MODEL. MASTER EQUATION. 

We assume in this paper the same atomic model as that 
presented in Refs. ||l^,|l^, i.e. a three-level A-system, 
composed of a ground-state level \g), a metastable state 
|e) and an auxiliary third fast-decaying state \r). Two 
lasers excite coherently the resonant Raman transition 
\g) —> |e) (with an associated effective Rabi frequency 
f2), while the repumping laser in or off-resonance with 
the transition \g) — ^ \r) pumps optically the atom into 
\g). With this three level scheme, one obtains an ef- 
fective two-level system with an effective spontaneous 
emission rate 7, which can be easily controlled by vary- 
ing the intensity or the detuning of the repumping laser 
pj|. In the following we follow the same notation as in 
Refs. |16|| . Let us introduce the annihilation and cre- 
ation operators of atoms in the ground (excited) state 



and in the trap level m [1), which we will call g™, 5^ 
(e/, e\). These operators fulfill the bosonic commutation 

relations [ffmisL] = ^rnm' and [eue],] = 6w . Using stan- 
dard methods of the theory of quantum stochastic pro- 
cesses pSJ-pql one can develop the quantum ME which 
describes the atom dynamics pi| 



p{t) = Cqp + Cip + C2P, 



(1) 



where 



Cap = -iHeffp{t) + ip{t)Hljj + Jp{t), (2a) 

Cip=-i[Hias,p{t)], (2b) 

C2P = -i[H,oh p{t)l (2c) 
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where 27 is the single-atom spontaneous emission rate, 
fl is the Rabi frequency associated with the atom tran- 
sition and the laser field, -qimikh) — (e, ^|e*'^^''|g, m) are 
the Franck-Condon factors, VV(0, 4>) is the fiuorescence 
dipole pattern, w^ (wf ) are the energies of the ground 
(excited) harmonic trap level m {I), and 5 is the laser 
detuning from the atomic transition. The new term re- 
spect to what is considered in Refs. [Q, is that of Hcoii, 
which describes the two-body interactions in the Bose 
gas. Only ground-ground collisions are considered be- 
cause we assume that the laser interaction is sufficiently 
weak to guarantee that only few atoms are excited (for- 
mally we consider only one). In the regime we want to 
study, only s-wave scattering is important, and then: 
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d^xip* ip* ip^n^-^^i, (4) 



where V'm denotes the harmonic oscillator wavefunctions 
and a denotes the scattering length. Following the sim- 
plifications of the QBME |l9|J2^ , we exclude the spatial 
dependence and therefore no transport or wave-packet 
spreading terms appear in (m). 



In the following we are going to work in the so-called 
weak-condensation regime, where no mean-field efects 
are considered. This means that we consider that the 
typical energy provided by the collisions is smaller than 
the oscillator energy. As shown in pO| , in typical experi- 
ments this condition requires that the condensate cannot 
contain more than 1000 particles. We shall work thus be- 
low such limit. We shall also consider that fl ^ ou. Also, 
due to the Festina-Lente requirements, ft is in general 
smaller than the typical collisional energy. Therefore we 
can formally establish the hierarchy Cq ^ C2 > J^i- 

Let us define a projector operator V: 



VX ^^\n){n\{n\X\n), 



(5) 



and its complement Q = 1 — V, and \n) = 
\No, Ni, . . .; g)'^\0, 0, . . . ; e) are the ground state config- 
urations with Nj atoms in the j'-th level, and no excited 
atoms. It is easy to prove that: 



CqV = VCiV = VC2V = 0. 

Projecting the ME (|l]), one obtains: 

v = V{Co + Ci+C2)w, 

w = QiCa +Ci+ C2)v + Q{Ci + £2)^, 



(6) 



(7a) 
(7b) 



where v ~ Vp and w — Qp. Laplace transforming 
{v{t) — > v{s)) and solving the system of equations (as- 
suming for simplicity w{Q) — 0), one obtains: 

sv{s)-v(Q)=V{Ca+Ci+C2) 

X [s - Q{Co + £1 + C2)]-\Ci + C2)i{s). (8) 

Performing the inverse Laplace transform, using (ph and 
V e-icp[— Cot]CiV = 7'exp[— £o''']'C2'P — 0, and applying 
the Markov approximation following Ref. [ll^, the ME 
becomes up to order 0{C\): 



V{t) = CcoUv{t) + Ccoolv(t) 



where 



Cc 



-VC2C0C2 



(9) 



(10) 



describes the collisional part. In principle in Ccoii appears 
a second term coming from the term between brackets 
in Eq. (pi), but since we consider the weak-condensation 
regime, we can employ the Born approximation as in Ref. 
p9[ to neglect these terms in the collisional part. There- 
fore the collisional part is described by a QBME as that 
of Refs. [^ 20|. The laser-cooling dynamics is described 
by 

+ -VLq \ dr \ dT'e^^^f^^^'^rie-^o^'ri, (11) 
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which has the same form of the ME calculated for the 
case of the laser cooling without collisions Q . 

Summaryzing, the dynamics of the system splits into 
two parts, (i) collisional part, described by a QBME, 
and (ii) laser-cooling part, described by the same ME as 
without collisions. The first correction to such splitting 
between both dynamics is of the order L2C\\ therefore 
the independence between the collisions and laser cooling 
dynamics is only valid in the weak-interaction regime. 

The independence of both dynamics, allows for an easy 
simulation of the laser cooling in presence of collisions. 
In particular, we simulate both dynamics using Monte 
Carlo methods, combining the numerical method of Ref. 
pot , with the simulations already presented in Refs. [ p^ . 
In the following two sections we shall present the results 
for one-dimensional and three-dimensional simulations 
respectively. 



III. ONE-DIMENSIONAL RESULTS 

This section analyzes the case of a one dimensional har- 
monic trap, and it is mainly devoted to the analysis of the 
laser cooling into states different that the ground state 
of the trap ||lj,|l^ (the ground-state case in analysed in 
the next section for the more interesting case of three 
dimensions). The laser-cooling into an excited state of 
the trap is only calculated in the one-dimensional case, 
because its analysis becomes very complicated in higher 
dimensions. The reason for that is that the ergodic ap- 
proximation, which we employ in Sec. IV, is incompatible 



with the analysis of the cooling in excited states of the 
trap. 

In this section we shall show that, as expected, the 
excited-state cooling is strongly affected by the collisions, 
even for modified low scattering lengths, in particular be- 
cause the collisions act as a mechanism to empty the de- 
sired excited state, and therefore compete with the cool- 
ing mechanism which tends to populate such state. For 
usual scattering legths and atom densities, the collisional 
processes are much faster than the typical cooling time, 
and therefore the excited state cooling is completely sup- 
pressed. However, the use of Feshbach resonances [ [l7[ , 
or lasers [18| , can modify the scattering length, in such a 
way that the typical time between collisions can be com- 
parable with the typical cooling time. In this section, 
we study the population dynamics for different modi- 
fied scattering lengths, analysing the transition from an 
ideal-gas regime (with a = as that studied in Refs. 
|16| ), to a regime in which atom-atom collisions consti- 
tute the dominant process. 



A. Collisional probabilities 

Following Ref. (2|] (Eq. (14)), the probability of a 
collision between two atoms respectivelly in the states n\ 



and 712 of an harmonic trap (of frequency w), to produce 
two atoms in the states n^ and 714 respectively, is given 
by: 

Pcou{ni,n2 -^ n3,n4) = 
— \U P 

■xNnANn^-5 



with 



J(iV„3+l)(7V„,+l + (5„3,„J, (12) 



where Nn denotes the occupation number of the level n 
of the trap. The Kronecker deltas in the previous ex- 
pression account for the bosonic factors appearing when 
two atoms in a particular level are created, or destruc- 
ted. Considering a highly anisotropic trap of frequencies 
u)x — iOy — Aw, Wz = oj, the problem becomes one- 
dimensional, and the expression (03) takes the form 

Pcou{ni,n2 -^ 713,724) = 

^wAc(7ll77, n2, 713)5(711 + 712 - 773 - TLi) 

xNn,{Nn, - Sn„n,){Nn, + l){Nn, + 1 + S„,.iJ, (13) 

where the Dirac S accounts for the energy conservation. 
In the previous expression: 



Ae(77l, 772, ^3) = [22("l+"^)77l !772!773! (77l + 772 - ^3)! 

2 
dzHn^(z)Hn^{z)Hn3{z)Hn^+n2-nAz)e^'^'' 



, (14) 



where Hj{z) is the Hermite polynomial of order j, and 



e 



16A2 / a 



n \aHo 



(15) 



with uho = \fhjrwjj. For typical experimental situa- 
tions a/oifo ^ 10"'^, so^ = 5xl0~6r^, where r = \a/ao, 
where Oq is the scattering length without any external 
modification. Therefore r = accounts for an ideal gas, 
whereas r — \ accounts for an unmodified scattering 
length. 



B. Laser— cooling transition probabilities 

The probability of laser-induced transition of an atom 
from the state ni to the state 772 of the trap is calculated 
in Refs. llq], and it is given by 



Pcool{ni -^ 712) ^ -—Ai{7li,7l2)Nni{Nn2 
27 



1) 



(16) 



where 



^2 r2TT i-TT 

Az(ni,772) = — / d^ desmew{e,(t>) 
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(17) 
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C. Dynamical cooling 



(18) 



Let us briefly review at this 
ing scheme proposed in refs. Hi 



oint the dynamical cool- 
|l6[ . It is easy to observe 
from the form of the rates (117), and from the same ar- 
guments as those used in Ref. [15|, that as in the single- 
atom case two different dark-state mechanisms can be 
employed: 

• "Franck-Condon" -dark-states. Let us assume a 
laser pulse with detuning S — slu respect to the 
atomic transition, where s is an integer number. 
It can be easily proved fl^ that a particular level 
of the trap |77i) remains unemptied (dark) if the 
Franck-Condon factor (777 -I- s| exp(7A:a;)|777) van- 
ishes. In particular, the dark-state condition for 
77 = 1 is 



77 = s 



1; 



(19) 



• "Interference"-dark-states. These dark-state 
mechanism is characteristic for dimensions higher 
than one. Let us assume the two-dimensional prob- 
lem, in which we have two orthogonal lasers char- 
acterised by two different Rabi frequencies: 51 in 
direction x, and Afl in direction y. The factor A in- 
dicates a possible difference between the intensities 
or phases of both lasers, and can be used to create 
a dark-state. If the laser detuning is zero and if we 



(777° le'''-" 1777: 



■S)/(^ 



^0\p^ky\ 
,0 



choose a value A = 

for a particular two-dimensional state 

then the selected level remains dark respect to the 

laser pulse 15 1. 






In absence of collisions, the above mentioned dark-state 
mechanisms [p5 16 1 allow to cool the atoms not only into 



the ground state, but also into an arbitrary excited state 
of the trap. In order to achieve that one has to use differ- 
ent dynamical cooling schemes. Each dynamical cooling 
cycle must contain sequences of pulses of appropriate fre- 
quencies. The following types of pulses are employed: i) 
confine7ne7it pulses: spontaneous emission may increase 
each of the quantum numbers TUx^y^z by 0(j]^). In D- 
dimensions pulses with detuning S — —Drfu)^ where if 
is the closest integer to rf ^ have thus an overall cooling 
effect, and confine the atoms in the energy band of D re- 
coils; ii) dark-state cooling pulses: these pulses should ful- 
fill dark state condition for a selected state to which the 
cooling should occur; iii) sideband and auxiliary cooling 



pulses: in general, dark state cooling pulses might lead to 
unexpected trapping in other levels. In order to avoid it, 
auxiliary pulses that empty undesired dark states and do 
not empty the desired dark state are needed; iv) pseudo- 
confining pulses: with the use of pulses i)-iii) cooling 
is typically very slow. In order to shorten cooling time 
we use pulses with 5 = —'irf'uj/2 and 5 — —rf'tu, which 
pseudo-confine the atoms below n — irf /2 and n = rf. 



D. Numerical results 

In the following we simulate the dynamics of the atomic 
population in the different trap levels using standard 
Monte Carlo methods. We consider the case of an har- 
monic trap with Lamb--Dicke parameter 7y = 3. We as- 
sume 7 = 0.04w, ri = O.OScj (consequent with the Festina 
Lente limit), and a number of atoms N = 133 (well in 
the weak-condensation regime). The calculations have 
been performed taking into account 40 trap levels. As an 
initial condition, we assume in all the following graphics 
a thermal distribution with mean (n) =6. In order to 
compare with the calculations without collisions, we anal- 
yse the same cooling scheme into the level n = 1 of the 
trap, as that studied in Refs. |1^], for the case without 
collisions. We consider cycles of four laser pulses with 
detunings 5 = slo, where si_2,3,4 = —9,8,-10,-3, and 
time duration T = 2j/fl^. Pulses 1 and 3 are confining 
pulses, pulse 2 is a dark-state pulse for ti = 1, and pulse 4 
is an auxiliary pulse. The one-atom emptying rates p5[ , 
|(n-f s| exp{ikx)\n)\'^ , for the first 10 levels of the trap are 
presented in Fig. n As one can observe, the effect of the 
pulses is to empty all the states except n — 1, which acts 
consequently as a trapping state. Observe, that due to 
the characteristics of the Franck-Condon factors, some 
levels of the trap are barely emptied, in particular for 
this case, n = 7 is also a quasi-dark state for pulse 2, 
and is poorly emptied by the auxiliary pulse 4. This 
is not important in the case without collisions, because 
n = 1 remains the darkest level throughout all the dy- 
namics. We shall show in the following that this is no 
more true when the coUions are accounted for. 
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FIG. 1. One-atom empty rates as a function of the trap 
level for the case of eta = 3.0, and a laser pulse with detuning 
S = SCO, with s = 8 (diamonds) and s — —3 (circles). 
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FIG. 2. Averaged population of the trap levels, for the 
case of r; = 3.0 and cooling cycles consisting of four pulses 
with durations T = 27/n^, and detunings 5 = slu with 
s = -9,8,-10,-3. The cases of (a) r=0.0, (b) r=0.4, (c) 
r=0.8, (d) r=1.2 and (e) r=5.0 are considered. 



Fig. H shows the evolution of the averaged population 
of the level n = 1 when r grows. We have evolved the 
system under 5000 cooling cycles (to avoid the effects of 
the initial conditions), and performed the average from 
the cycle 5000 until the cycle 15000. In Fig. we have de- 



picted the averaged population distribution for the cases 
of (a) r = (ideal gas), (b) r = 0.4, (c) r = 0.8, (d) 
r = 1.2, and (e) r = 5.0. For the ideal-gas case r = 
one obtains that the population is completely condensed 
into the level n = 1 |16|. When r is increased the laser- 
induced condensation into the excited state n = 1 is de- 
stroyed, but in a non-trivial way. We observe in Fig. 
(b) that for r = 0.4 the population is basically distributed 
in two well defined peaks, one in n = 1 and the other in 
n = 7. The reason for this behavior can be understood 
very well using Fig. Il|. In absence of collisions the level 
n = 1 is not emptied at all, while, as pointed out pre- 
viously, n = 7 is emptied, but slowly, and therefore at 
the end the population is finally transferred to the level 
n = 1. However in the presence of collisions, level n — 1 
is still a dark-state for the laser, but it is emptied by 
the collisions with a frequency proportional to Ni. This 
means that the population is pumped into n — I due to 
the laser cooling, but the more the population we pump 
into n — I the more the level is emptied via collisions. 
This effect can be well illustrated by Fig. ||, where one can 
observe periods of filling of n = 1 followed by abrupt de- 
cays of the level population. The emptying of level n = 1 
is mainly produced via collisions between two atoms in 
the level n = 1 to produce two atoms in n = and n — 2 
respectively. The laser cooling provides a mechanism to 
repump such expelled population from the level n — 
and n — 2 back to the level n — 1. Such control is al- 
ready maintained for large occupations of n = 1, but in 
an unstable way, due to the highly non-linear character 
of the dynamics. A slight excess of population into the 
level n = and n = 2 provoques a speed-up of the emp- 
tying process of ri = 1. This situation is reflected, for 
example, in the behavior of the system between cycles 
4500 and 5500 in Fig. pi In particular, level n = \ can 
become more emptied than n = 7, and the latter turns 
to be the effective darkest level, i.e. the level less emp- 
tied. Therefore the population tends to be transferred 
into n = 7. But, when Ni increases so does the empty 
rate of the level n — 7, which can become larger than 
that of n = 1, and so on. Therefore, as consequence of 
this process a non-linear pseudo-oscillatory motion be- 
tween the populations of n = 1 and n = 7 is produced, 
as observed in Fig. pi This oscillatory motion leads to 
the two-peaked distribution of Figs. 0. Finally, when r 
becomes very large the collision dynamics is much faster 
than the cooling time, and the peaked structure dissa- 
pears, as observed in Figs. || (c), (d) and (e). Observe 
that nevertheless the efl'ects of the laser cooling mecha- 
nism are nevertheless present in Fig. ^ (f). In absence of 
laser cooling, it can be demonstrated that for r — 5, the 
population has a maximum in n = 0. On the contrary, 
in the presence of the laser, the population of n = is 
very efficiently and rapidly emptied by the pulse with 
detuning s = 8, which is the most rapid cooling process. 
For larger r even this process is eventually overcome, and 
one recovers the same distribution as that obtained only 
considering the collisions without laser cooling. 



IV. THREE-DIMENSIONAL RESULTS 

In this section we analyze the case of the laser-cooling 
into the ground state of an isotropic three-dimensional 
harmonic trap, of frequency lu. The numerical calcula- 
tion of the system dynamics for the three-dimensional 
case is quite complicated, due to both the degeneracy of 
the levels, and the difficulties to obtain reliable values 
for the integrals C/(ni, 712, 713, 714). Therefore, we shall 
limit ourselves to the use of the ergodic approximation, 
i.e. we shall assume that states with the same energy are 
equally populated. The populations of the degenerate en- 
ergy levels equalize on a time scale much faster than the 
collisions between levels of different energies, and than 
the laser-cooling typical time. This approximation leads 
to the correct steady-state distribution, although the dy- 
namics can be slightly different than in the non-ergodic 
calculation EOfl. 
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FIG. 3. Detail of the dynamics of the populations of lev- 
els n = 1 (solid line) and n — 7 (dashed line), for the case 
of 77 = 3.0, r — 0.4. Each cooling cycle consists of four 
pulses with durations T = 2^/QF', and detunings 5 — suj 
with s = —9,8,-10,-3. The non-linear pseudo-oscillation 
between n = 1 and n = 7 can be clearly observed. 



Following ref. ||29[] the probability of a collision of two 
atoms in energy shells rii and n2, to give two atoms in 
shells 713 and 714 (where this collision is assumed to change 
the energy distribution function) , is of the form: 

P(7ii, 7^2 -> 7i3, 7^4) = A(7ij -I- l){nj + 2) 

X . 



(20) 



where g„^ = {nk + l)(nfc + 2)/2 is the degeneracy of 
the energy shell ti^, rij = 777,m{7ii, 712, 713, 714}, and A — 
{Aa?'Ld'^ra)/{TTh) ~ 1.5 x 10^^. Concerning the laser- 
cooling probabilities we shall use the same expressions 
as those already developed in Refs. |16|] . 

In the following we simulate the evolution of the sys- 
tem by using again Monte Carlo simulations. Due to 
numerical limitations we consider a Lamb-Dicke param- 
eter 77 = 2. We assume as previously 7 = 0.04a; and 
il = O.OSw, and a number of atoms N = 133. As a 
first step, we begin with a thermal distribution of mean 
{n) — 6, and evolve the system just with collisions, un- 
til obtaining a Bose-Eintein distribution (BED) (which 



does not coincide exactly with the thermodynamical one, 
due to finite-size effects), see Fig. ^. The distribution 
obtained in this initial step serves as the initial state for 
laser cooling. As we see, it already contains quite subtan- 
tial amount of atoms condensed in the ground state, but 
also a lot of uncondensed ones. Laser cooling will trans- 
fer the latter ones into the ground state. We apply our 
laser cooling cycles, each one of them composed by two 
laser pulses of detuning S = scu, with si,2 — —4, 0, and 
time duration T = {2^)/Vt^ . The laser pulses are emitted 
in three orthogonal directions x, y and z, and are charac- 
terized by their respective rabi frequencies ^^ = ^y = ^, 
fiz = Az^. For the first pulse we assume Az — 1, while 
for the second one Az = —2 is considered. With this 
choice, the second pulse is an "interference" -dark-state 
pulse for the ground-state of the trap. Fig. @(a) shows 
(dashed line) that these two pulses are able to condense 
the population into the ground state of the trap, in ab- 
sence of collisions; in particular no confinement pulses (of 
detunings S = —12uj in this case) are needed. This is due 
to the bosonic enhancement and the fact that initially 
the system is already partially condensed. The dark- 
state pulse is neccesary to repump the population in those 
states of the energy shells 1, 2 and 3, which are dark re- 
spect to the pulses with detuning S = —Alu. Fig. || shows 
(solid line) the dynamics of the population of the ground- 
state in presence of collisions. After 600 cycles, all the 
population is transferred to the ground state of the trap. 
This means that applying the laser cooling scheme brigs 
the system into an effective BED of T = 0. It is easy to 
undertand why the effect is maintained in presence of col- 
lisions, even considering that the coUisional dynamics is 
much faster than the laser-cooling one. The laser-cooling 
mechanism tends to decrease the energy per particle (i.e. 
the chemical potential of the system) , in the same way as 
evaporative cooling does, but without the losses of par- 
ticles in the trap during the process. Thermalization via 
collisions brings the system to a lower temperature. Re- 
peating the laser cooling sufficient times the system ends 
with an effective zero temperature. Finally, let us point 
out that some auxiliary pulses which are needed in the 
ideal gas, are not in presence of collisions. In particular 
for the previous example, the pulse of zero detuning (re- 
quired for the ideal gas case, FigM5|(b) dashed line) is no 
more needed, as shown in Fig. pTb) (solid line). Thus, 
the laser-cooling scheme is not only possible in presence 
of collisions, but can be even significanly simplified. 
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FIG. 4. Average of the population of the different energy 
shells, after evolving the system under collisions, starting from 
a themal distribution of (n) = 6. 
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FIG. 5. Dynamics of the population of the ground state of 
the trap for the case of ry = 2, 7 = 0.04u>, Q = O.OScj, N = 133, 
starting from the BED distribution of Fig. kl (a) Cooling 
cycles composed of two laser pulses of detunings S — suj, with 
Si,2 = —4,0, and A^, = 1,-2 repectively, are used. The 
simulations taking into account the collisions (solid line) and 
no taking them into account (dashed hne) are compared, (b) 
Same case, but only employing pulse 1. 



V. CONCLUSIONS 



In this paper, we have analysed the effects of the atom- 
atom collisions on the colective laser cooling of bosonic 



gases trapped in an harmonic trap, under the Festina- 
Lente condition. In particular, we have studied the case 
in which the mean-field energy provided by the atom- 
atom collisions is much smaller than the typical energy 
of the harmonic trap. Under such conditions, we have 
derived the ME which describes the system, and ob- 
served that such ME splits into two parts: (i) a purely 
collisional part which has the form of a QBME, and (ii) 
a purely laser-cooling part, which has the same form as 
the ME which describes the laser-cooling in absence of 
collisions. By using this ME, we have simulated the dy- 
namics of the trapped gas for different situations. First, 
we have analysed the cooling into an excited state of a 
one-dimensional trap. We have observed that the tran- 
sition from the ideal-gas limit (in which the atoms are 
completely condensed into the chosen excited state) to 
the case in which the collisions dominate the dynamics, 
is not trivial, specially when the collisional and cooling 
time scale are comparable. In such a case, cooling and 
collisional processes enter in competition, and new phe- 
nomena can appear, as for example unstable population 
of an excited state followed by abrupt population decays, 
and non-linear pseudo-oscillations between different trap 
levels. We have finally analyzed the laser-cooling into 
the ground state of an isotropic three-dimensional har- 
monic trap, by using the ergodic approximation. We 
have shown that, although the collisional time is typically 
much faster than the cooling time scale, the laser cooling 
allows to transform a BED with a finite temperature into 
an effective BED with zero temperature. The laser cool- 
ing reduces the chemical potential of the trapped atoms, 
while the collisions provide the thermalization. 

Let us finally present important remarks concerning 
the scaling of our theory, the situation beyond the weak- 
condensation regime, and the problem of the two- and 
three-body losses in the trap. First, we stress that we 
have presented here the results obtained for 77 = 2, 3 
only for the reasons of numerical complexity which grows 
rapidly with rj. Qualitatively, the same results can be ob- 
tained for larger ry's, and therefore, for lower densities. In 
fact, we have observed similar results for 77 = 5 in one di- 
mensional simulations. If we increase r] by factor F, the 
corresponding density (for fixed A'^) decreases as F~^, the 
three body loss rates as F~^, whereas the trap frequency 
decreases as -F^^, which means that the corresponding 
cooling time (to fulfill the Festina Lente conditions) will 
increase as F^, i.e. much less than the lifetime due to 
three-body collisions. 

Beyond the weak condensation regime, the mean-field 
energy cannot be neglected, and therefore the trap levels 
are no more the harmonic ones. This has a two-fold con- 
sequence: (i) The levels of the trap are non-harmonic, 
i.e. they are not equally separated, because their ener- 
gies become dependent on the occupation numbers; (ii) 
the wavefunctions are different, and in particular the con- 
densate wavefunction becomes broader (we consider here 
only the case of repulsive interactions, a > 0). The fact 
that the energy levels are not harmonic any more, com- 



plicates the laser cooling, but the use of pulses with a 
variable frequency and band-width should produce the 
same results as those presented here. The point (ii) im- 
plies that the central density of the interacting gas is 
much lower than the one predicted for noninteracting 
particles. In fact, the ratio between the interacting-gas 
central density (in Thomas-Fermi (TF) approximation) 
and the ideal-gas central density, goes as ISOl 



riTF 152/57ri/2 / Na 



'^ideal 



a-HO 



-3/5 



(21) 



where the central density for the ideal case is given by 
n.Meai = A^/(7r3/2a|^o). 

The above result has important consequences, when 
one considers the problem of three-body collisions, which 
usually begin to play a role at densities of the order of 
10^^ atoms/cm'^. For example, let us analyse the case 
of Sodium, for which an = ^'h/raiOR = 0.132/im, and 
a = 2.75nm. From the definition of 77, ano = V'^r- For 
the ideal gas case, the regime in which three-body losses 
are important is reached for N ~ 12.877'^. For 77 = 8 
this means N — 6.5 x 10'^. Amazingly, for the interact- 
ing gas, the same is true for N ~ 5.177^, and therefore 
for rj — %^ the regime in which three-body losses are 
important is reached for A^ = 1.3 x 10^. Below this num- 
ber, the interaction between the particles is dominated 
by the ellastic two-body collisions considered in this pa- 
per. As point out above, our laser cooling scheme could 
be extended beyond the weak-condensation regime, and 
therefore laser-induced condensations of more than 10^ 
atoms are feasible. 

Concerning other loss mechanisms, we have to mention 
here the hyperfine changing two-body collisions, or gen- 
erally speaking any inelastic two-body processes. These 
can be supressed completely if we cool atoms into the 
absolute ground internal state, which is possible in the 
dipole traps. For alkahs this is typically done by coohng 
to the lowest energy state in the lower hyperfine mani- 
fold (external static magnetic fields are used to split the 
levels within the hyperfine manifold). 

Finally, yet another loss mechanism disregarded here 
is due to photoassociation, i.e. excitation of molecular 
resonances. This kind of loss rates are typically of the 
order ^n{X/2'KY (where 7 is the linewidth of the aux- 
iliary level I?'), A is the laser wavelength, and n is the 
atomic density), i.e. allow for achieving about 1000 cool- 
ing cycles of duration ~ I/7 provided the density remains 
smaller than 10^^ — 10^^ atoms/cm^. However, the pho- 
toassociation losses can be reduced by several orders of 
magnitude if the laser is red detuned, and tuned exactly 
in the middle of the molecular resonances ||3^] (note that 
this is the detuning respect to the one-photon transitions 
from the ground states to the state |r), and not the two- 
photon detuning). Other, possibility, is of course to use 
a more intense laser tuned below the Condon point, i.e. 
the minimum of the molecular potential. 
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